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In this talk we describe our approach for the computation of multi-leg one-loop amplitudes and present some 
first results relevant for LHC phenomenology. 



1. Motivation 



In the next year the Large Hadron Collider at 
CERN will start operating at the up to now un- 
reached TeV energy scale. Apart from putting 
our theoretical understanding of the electroweak 
symmetry breaking mechanism under scrutiny, 
new effects beyond the Standard Model that oc- 
cur in this energy region can also be explored 
|1I2| . A successful quantitative analysis and com- 
parison of the LHC data with theoretical predic- 
tions will crucially rely not only on precise predic- 
tions for a variety of signal cross sections, but also 
on a sufficiently precise knowledge of the various 
backgrounds. Although in some favourable cases 
the Standard Model backgrounds can be mea- 
sured directly, for other important search chan- 
nels this is not possible and backgrounds have 
to be determined through theoretical predictions. 
This is especially worrisome, as QCD predic- 
tions for multi-particle processes are plagued by 
large uncertainties arising from high powers of as . 
Large scale uncertainties make it mandatory to 
include radiative corrections beyond the leading 
order in QCD. Unfortunately, due to the com- 
plicated singularity structure of QCD and the 
general combinatorial complexity of multi-scale 
problems, one-loop computations with more than 
four external partons/particles are a highly non- 



trivial and time-consuming task. Fortunately, at 
the dawn of the LHC era, this research direc- 
tion is very actively pursued by multiple groups 

|:mi6lVlxl9ll()llill:jli3ll4| . 

In this talk I present our approach to this prob- 
lem. The aim of our collaboration is to build 
a flexible and reliable tool, which allows for the 
evaluation of one-loop multi-leg amplitudes. It 
is based on a combination of algebraic and nu- 
merical methods. The corresponding project is 
called GOLEM (General One-Loop Evaluator for 
Matrix elements) and will be discussed in the 
next section. Section 3 is dedicated to first appli- 
cations of these methods relevant for LHC phe- 
nomenology. 

2. The GOLEM project 

A General One-Loop Evaluator for Matrix ele- 
ments (GOLEM) should proceed from a Feynman 
diagrammatic representation of a given scattering 
amplitude to a computer code which provides a 
numerically stable and accurate answer for the 
desired cross section. 

In a first step each amplitude can be repre- 
sented in terms of irreducible Lorentz tensors 
which correspond uniquely to gauge invariant op- 
erators: 
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The mapping of the diagrammatic input onto 
such a tensorial basis can be accomphshed with 
algebraic manipulation programs. We generally 
use FORM 3.1 ^^l- Each Feynman diagram Qg 
can be projected onto these sub-amplitudes: 
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The coefficients C of the gauge invariant operators 
are momentum space integrals where free Lorentz 
indices are contracted with external momentum 
vectors. In a third step these coefficients, which 
depend on Mandelstam variables Sij — (pi +Pj)'^ 
only, are expressed in terms of some integral basis 
to be discussed below: 
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In a forth step, the coefficients have to be ex- 
ported to some numerical programming language. 
Optionally they can be simplified beforehand us- 
ing algebraic programs like MAPLE or MATHE- 
MATICA: 



A{\pj),, 



BIG 

X simplify[CB/G(sjfe, 
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We now discuss the integral basis used for our 
reduction algorithms. Following '16' one can rep- 
resent any given tensor integral by Feynman pa- 
rameter integrals in shifted dimensions: 



By applying differentiation by parts in parameter 
space and d = 4 kinematical identities ^ one 
can map all these integrals to the following basis 
without the need for higher dimensional scalar 
integrals with iV > 4: 
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and scalar integrals /g", /;^'+^ /^'+^ 
call this basis the GOLEM basis. Note that 
any iV-point function can be mapped to this 
basis without introducing inverse Gram deter- 
minants, i.e. l/det(G)'' with Gij — 2ri ■ rj, 
Tj — pi + . . . -\- pj^ which occur in Passarino- 
Veltman reduction. There are three alternatives 
for evaluating the basis elements: 

1. algebraic reduction to the one- loop master 
integrals /4+^ 

2. semi- numerical reduction to scalar 
integrals, 

3. direct numerical evaluation. 

In this context, semi- numerical means that the 
same reduction formulas as in 1. are applied at 
the numerical level. In case 1. and 2. Gram de- 
terminants are reintroduced. The direct numer- 
ical evaluation is chosen whenever numerically 
critical phase space regions are approached. 

For the numerical evaluation of the GOLEM 
basis integrals two different methods have been 
developed. The first is based on contour defor- 
mation in Feynman parameter space 7 and the 
second on performing the Cauchy integration ex- 
plicitly The latter method is faster, since 
(at most) 2- rather than 3-dimensional represen- 
tations are used. We have tested the robustness 
of the numerical routines for large samples of re- 
alistic phase space points. 

The algebraic/numerical algorithms are imple- 
mented in the flexible Fortran 90 code G0LEM90. 
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3. Applications to LHC phenomenology 

Using the algorithms discussed above, several 
computations relevant for LHC phenomenology 
have been carried out so far. We present three 
examples here, others can be found in [17118119] . 

3.1. The process gg W*W* IvVv' 

This process contributes to the dominant irre- 
ducible background to the Higgs production pro- 
cess gg ^ H ^ W*W* . It is also relevant below 
the WW threshold. The gluon-vector boson in- 
teraction is mediated by a quark loop. On-shell 
production results with massless ^ and massive 
quarks |21| have been known for a long time. For 
on-shell W bosons, a phenomenological analysis 
can be found in |22j. We have now a full calcu- 
lation for off-shell VF-pairs, where the first two 
quark generations are treated massless, but the 
bottom and top masses in the third generation 
are kept. In the latter case the amplitude con- 
tains six different scales: the Mandelstam vari- 
ables s, t, the ly-virtualities S3, S4 and the quark 
masses m;,, mf. For this amplitude, it is pos- 
sible to reduce algebraically down to master in- 
tegrals, J^, /g , I"^^. The amplitude is repre- 
sented by nine independent gauge- invariant struc- 
tures. Each coefficient can be written as a lin- 
ear combination of 27 basis functions. The al- 
gebraic simplification results in amplitude repre- 
sentations which have at most one inverse Gram 
determinant. The Gram determinant is related 
to the transverse momentum of the W bosons 
by p^{yV) = det (G)/s^. In this computation 
there is a tiny phase space region pt{W) < 0.1 
GeV, 1 53^4 — M^l Mw^w which is numeri- 
cally problematic, but it does not contribute siz- 
ably to the relevant cross sections 23 . In Tab. ^ 
we show typical LHC cross sections for two mass- 
less generations and including the third genera- 
tion. We also compare to the full NLO result for 
qq —> W*W*, which is computed using MCFM 
|24|. Cross sections without selection cuts (tot), 
with standard LHC cuts (std) and Higgs search 
cuts (bkg) are given (see for details). This 
calculation shows that for Higgs search cuts the 
gg W*W* process enhances the previously 
known T4^-pair background by approximately 30% 





gg (2 gen.) 


gg (3 gen.) 


<7N LO+<7gg 

aNi.n 




53.61(2) 


60.12(7) 


1.04 




25.89(1) 


29.79(2) 


1.06 


CTbkg 


1.385(1) 


1.416(3) 


1.30 



Table 1 



Background cross sections in fb for charged 
lepton pair production in giuon fusion gg —^ 
W*W* Ivl'v'. Results for two and three 
quark generations are shown for different selec- 
tion cuts. The last column shows the importance 
of the gluon fusion process (3 gen.) relative to 
the quark induced channel in NLO QCD. 



and has to be taken into account in experimental 
studies. In Fig. we show the charged-lepton 
azimuthal opening angle distribution. 

3.2. The process gg HHH 

Multi-Higgs production via gluon fusion allows 
in principle to measure the triple and quartic 
Higgs self-coupling. A full calculation of three 
Higgs boson production in gluon fusion was ac- 
complished recently by Four topologies 
with different sensitivities can be distinguished. 




The amplitude can be decomposed into five in- 
dependent gauge invariant structures iT{TiT2)^ 
P2- ^1- Pi Pi- ^2- Pj, I, j S {3, 4}, where is the 
gluon field strength tensor. Some of the respec- 
tive coefficients are related by Bose symmetry. 
Evaluating all coefficients and testing these re- 
lations serves as a check of the calculation. Alge- 
braic reduction to 68 master integrals is possible 
and one finds a numerically stable representation 
of the amplitude with at most one inverse Gram 
determinant. The result is summarised in Fig. 
121 Destructive interference effects are present be- 
tween pentagon box and box triangle topologies. 
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Figure 1. Azimuthal opening angle distribution 
for the charged lepton pair at the LHC. The upper 
(lower) solid line is the gluonic contribution (x 10) 
including (excluding) the massive third family 
quarks. The other curves are the Higgs signal 
(x5) and LO/NLO quark contribution (MCFM). 
For details see |23) . 



The total cross section for mn = 120 GeV, being 
only of the order of 0.06 fb, is out of reach at the 
LHC. Note that in Beyond Standard Model sce- 
narios with an enlarged fermion or Higgs sector 
this cross section could be much larger. We com- 
pared our full result to the heavy top limit and 
observed that for three Higgs boson production 
the mtop — > oo approximation is not applicable. 
Here, the kinematically favoured region is above 
the 2 rritop threshold, rather than below, as re- 
quired by the heavy top limit. 

3.3. The PP 4 jet amplitude 

Up to now no full LHC process with 2-^4 
kinematics has been evaluated. To have a good 
theoretical prediction for the huge four jet rate, 
next-to-leading order corrections to the partonic 
reactions gg gggg, gg ggqq, gg qqqq, 
qq — > qqqq plus crossings need to be evalu- 
ated. By using our GOLEM algorithm we have 
produced a representation of the six quark am- 
plitude in terms of the GOLEM basis integrals 
([T|) which allows for the semi-numerical evalu- 
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Figure 2. The total triple Higgs production cross 
section vs. mn for different topologies (tril ^ 
X4H, tri2 Ag^). The H — helicity component 
and all triangle topologies are suppressed. 



ation of this amplitude. The evaluation time 
of one kinematical point is of the order of a 
second. Alternatively we have produced an al- 
gebraic representation of the amplitude in terms 
of master integrals /g and I^'^'^. The two 
independent approaches allow for an efficient de- 
bugging of our algebraic/numerical procedures. 
As this work is not yet completed we just illus- 
trate this by providing a numerical evaluation 
of the most complicated Feynman diagram in 
this calculation, the 6-quark pentagon diagram 




This diagram can be written as 
A^i-^^^^-^''^=^«(fci,...,fc6) = 

For the helicity amplitude and the 

kinematical point 
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7000. 
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0. 


7000. 
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7000. 


0. 


0. 


-7000. 


ks 


5105.37 


1752.05 


-972.19 


-4695.74 


ki 


4582.45 


-3465.91 


1373.64 


2664.47 


k5 


2917.90 


2182.21 


876.20 


1727.53 


ke 


1394.27 


-468.35 


-1277.65 


303.74 



we find, up to an overall phase, 

A = 3.27888-1- i 0.950546 , 
B = 10.7288+ a5. 53100 , 
C = -1.16693 i 59.1171 . 

The evaluation time for this point with our code 
is about 0.2 sec on a standard PC (Pentium 4, 
2.8 GHz). The evaluation of the full 6-quark am- 
plitude including real emission corrections is in 
progress. 

4. Conclusion 

In this talk we presented our GOLEM approach 
for the evaluation of one-loop multi-leg processes 
relevant for precise LHC phenomenology. The 
formalism allows in principle for the evaluation 
of general iV-point processes with massive and/or 
massless particles. The algebraic/numerical re- 
duction and evaluation techniques are imple- 
mented in the Fortran 90 code GDLEM90. As appli- 
cations, recent results for the gluon-induced pro- 
cesses gg W*W* IvVv', gg ^ HHH and 
the hexagon process qq — > qqqq were presented. 
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